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Abstract 

The correlations between the segments of a semidilute polymer solution are found to induce 
correlations in the positions of small particles added to the solution. Small means a diameter 
much less than the polymer's correlation length. In the presence of polymer the particles behave 
as if they attracted each other. It is shown how the polymer's correlation length may be 
determined from a scattering experiment performed on the spheres. 

1 Introduction 

Small colloidal particles in a semidilute polymer solution are shown to interact via a long ranged, 
attractive potential of mean force. This potential of mean force Q is an effective interaction 
potential between the particles due to the presence of the polymer. The particles are small with 
respect to the polymer's correlation length and, as it is this length which determines the range of 
the attractive potential of mean force, the attractive interactions are long ranged with respect to 
the direct interaction between the spheres. The effective attractions between the colloidal particles 
are a result of the correlations between the segments in a semidilute polymer solution Q. The 
segments of a polymer molecule in a good solvent follow a self-avoiding walk. Because of this the 
fraction of the volume of the solution occupied by the polymer segments is highly non-random: if 
there is a segment at one point then the average density of polymer nearby is much higher than the 
average density of segments. This implies that the volume not occupied by the polymer segments 
is also highly correlated; if you imagine a picture in which every point is either black or white, 
then if there is a pattern to the white areas then there must be a pattern to the black areas. But 
the particles are only free to move in this volume which is free of polymer; they are excluded from 
the volume occupied by the polymer. So the particles move in a fraction of the solution's volume 
which has long range correlations and so their positions have long range correlations. In fact the 
correlations between the spheres resemble those produced by a long ranged attraction. As the 
correlations in the spheres reflect those in the polymer, scattering experiments |^ on the colloidal 
particles may be used to probe the structure of a semidilute polymer solution; in particular they 
provide a way of measuring the correlation length of the polymer. 



^Present address: Department of Chemistry and Biochemistry, University of California 
Los Angeles, Los Angeles, CA 90024, USA. email: sear@chem.ucla.edu 
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The colloidal particles could be protein molecules [Q, surfactant micelles |^], or synthetic, 
polymer or silica, spheres HI, ^. Both protein molecules and micelles typically have diameters of 
a few nm. The direct interaction between particles is taken to be purely repulsive, there are no 
attractive interactions. The polymer is considered to be in a good solvent and is therefore swollen 
due to self-interactions The polymer's radius of gyration is much larger than the diameter of the 
particles, ~ lOnm or more. See Refs. ^ for examples of work where the radius of gyration 

is comparable to or smaller than than the diameter of the particles. Recent work of Khalatur et 
al. |13| has dealt with similar mixtures to those considered here using RISM integral equations. It 
is more numerical and less intuitive than here. The RISM approach is also best suited to higher 
polymer densities than those considered here, say volume fractions of polymer of 10% and more. 
Our simple scaling approach breaks down at these densities. All interactions are assumed to be 
excluded volume interactions, that is two polymer segments may not occupy the same volume and 
likewise for two particles or a particle and a polymer segment. This is reasonable if the solvent is 
good for the polymer, and the polymer does not absorb onto the particle. 

All three interactions are excluded volume and therefore there are no energy scales, apart from 
the temperature T. Therefore, the behaviour is a function only of the length scales of the mixture. 
The length scale of the particles is just their diameter D. A pure semidilute polymer solution has 
only one relevant length scale, the correlation length £^ |^, which is roughly the distance between 
interactions between segments on different polymer chains. The radius of gyration is not a relevant 
length scale as it is much larger than ^ and the chain loses its correlations and 'forgets' which 
polymer it belongs to over a distance of ^. The mixture then has only two relevant length scales: 
D and [14, 15, This means that the particle-polymer interaction depends only on their ratio: 



2 Theory 

The polymer molecules consist of flexible linear chains of segments of size a. These segments are 
much smaller than the spherical particles, D ^ a and so the volumes from which the segments 
exclude spheres overlap. Each segment excludes a sphere from a volume ~ (D + a)^ centred on itself 
but as the segments are part of a chain there are always two segments a away from any segments, 
another two ~ 2a away, etc.. This problem of the excluded volumes overlapping can be reduced 
if the polymer segment length is rescaled from a to D. This is illustrated in Fig. 1, which shows 
a sphere and a part of a polymer chain. As semidilute polymer solutions are scale invariant at 
length scales much less than the correlation length we are free to perform this rescaling. Now 
that the polymer segment length is D the excluded volumes of adjacent segments overlap much less 
strongly, although they still overlap. Each rescaled segment consists of many segments, for D ^ a, 
and so acts like a polymer coil of the same size as the particles. The particles cannot penetrate 
coils of their own size [17[ and so the segments exclude particles from a volume ~ D^. The ratio 



of the number of segments of size D to those of size a is ~ (a/D)^^^ if the polymer is in a good 
solvent. 

From now on we only consider the scaled polymer segments of size D, see Fig. 1. Treating the 
polymer as a sequence of rescaled segments of size D makes it clear that the only relevant length 
scales are D and the polymer's correlation length. So, the polymer density cp is the density of 
these polymer segments. These segments exclude each other from a volume ~ just as dilute 
polymer coils exclude each other from a volume of order of their radius of gyration cubed. We 
have approximated the interaction between a segment and a sphere by a hard-sphere interaction 
of diameter D; the same interaction as between the spheres. This is an approximation of course; 
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the segments do not interact with a completely hard potential, they are only very approximately 
spherical and they are connected into chains. 



2.1 Scaling theory for semidilute polymers 

First, we briefly summarise the scaling theory for polymer solutions in good solvents, it is described 
in detail by de Gennes |2[. Our spheres will be added to a semidilute solution of polymer molecules 
in a good solvent. As noted above we have rescaled the polymer segment length from a to D. The 
size of an single isolated polymer coil oi ud segments is measured by its radius of gyration Rq 
which is given by Rg ~ Dr?^ ^ in a good solvent Q. The correlations (fluctuations) in an isolated 
polymer coil extend across the entire coil; the correlation length ^ is thus ~ Rq. When polymer 
coils overlap, interactions between different chains reduce the correlation length ^. The overlap 
volume fraction c* is Cp = D~^nj^^^. This volume fraction is the boundary between the dilute and 
semidilute regimes; below c* the interactions between chains can be treated as a perturbation but 
above c*p the interaction between a pair of chains is ^ ksT. In the dilute regime ^ = Rg, and in 

the semidilute regime. = D'^^'^Cp^^ . 

We require the pair distribution function for the polymer segments gp{r), where r is the mag- 
nitude of the separation of the two segments [18|. It is. 



gp[r) ~ < 



r<D 

4/3 

D<r<C . (1) 

1 I gexp(-r/g) f < r 



Note that we define our pair distribution function differently from Ref. Q. The first line of Eq. 
(|^) just refiects the mutual impenetrability of polymer segments, the second is the density of a self- 
avoiding random walk, and the third line reflects the screening of correlations at distances greater 
than ^. Note that for r ^ ^ the dominant contribution to gp is the contribution from a single 
chain. 

2.2 Scaling theory for spheres in a semidilute polymer solution 

In order to examine the pair correlations and hence the scattering from spheres in a semidilute 
solution we require the pair distribution function of the spheres, gsir). For r < D this is due to 
the hard-sphere interaction but for r > D gs is not equal to 1 due to the fact that the cavities in 
the polymer solution are correlated; see The cavities are the points in the polymer solution 
where there is no polymer within a distance D, and so where it is possible to insert a sphere [ p!s[| . 
In the limit of very low sphere density cs we need only consider the two spheres which form the 
pair; there are no other spheres present in this limit and so the correlations between the spheres 
are entirely due to correlations in the cavities: 

lim gs{r) = | ° , , ^ n • (2) 

cs^o ^ ' \ ypir) r > D ^ ^ 



where yp{r) is the cavity-cavity correlation function equal to the probability of inserting a 
pair of spheres a distance r apart divided by the square of the probability of inserting one sphere. 
The probability of inserting a single sphere will be denoted by a. Below we start by estimating a 
and then go on to consider the probability of being able to insert a pair of spheres which yields an 
estimate oi yp. 
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The concentration of segments is cp and each segment excludes spheres from a volume of ~ D^. 
If we ignore the fact the segments are connected then as the polymer density is low, we can also 
neglect the overlaps of excluded volumes of nonadjacent segments and we arrive at the conclusion 
that the fraction of the volume denied to the spheres is cpD^ |14, 15|. Using ^ = D'^/'^Cp^^^ this 



fraction can also be written as (D/^)^/^, and the fraction of the volume available to the spheres, 
a, is 

a ~ 1 - - . (3) 



I) 

Equation (^) has a probabilistic interpretation which will be used below. It is the probability that 
there is no polymer segment within D of a randomly chosen point in the solution. 

In order to calculate the pair distribution function between spheres, gs, we relate it, using Eq. 
to the distribution function between pairs of cavities yp, i.e., pairs of points with no polymer 
segments within a distance D of either of them. Then we express yp in terms of a and gp, both 
of which are known, Eqs. (^) and @). In order to relate yp to gp we use the fact that for a pair 
of points, labeled 1 and 2, there are 4 possibilities: 1 and 2 both have polymer segments within 
a distance D; 1 does but 2 does not; 2 does but 1 does not; and neither does. Therefore, the 
probability of the fourth possibility occurring, which is a'^yp, is equal to 1 minus the sum of the 
probabilities of the other possibilities. The probability of both points being within D of a, polymer 
segment is the square of the probability of one point having a nearby segment, {D/(^)^^^, times 
gp{r), where r is the magnitude of the separation r between points 1 and 2. In deriving this result 
we have integrated over spheres of diameter D around both points; therefore, it is only valid for 
r > 2D but we will use it for r > D. This is a very minor approximation, yp still has the correct 
behaviour near r = D. The probability of point 1 being near a segment and point 2 not being 
near a segment, is the probability of point 1 being near a segment, (D/^)^/^, times the probability 
of there being no polymer segment near a point r away from a seg ment 1 - {D / S)^/^ gp{r). This 
second probability is just 1 minus the probability of there being a segment at a separation r from 
another segment. So, 

(l - {D/O^I^)\p{r) ~ 1 - {D/0^l^gp{r) - 2{D/0'''' (l - {D/ifl^gp{T)) r > D, (4) 
and 

Note that a^yp{D) ~ a which is at it should be, the probability of inserting two particles next to 
each other is almost the same as the probability of being able to insert one particle. The probability 
of inserting a sphere close to where a sphere has already been successfully inserted is higher than 
at a randomly chosen point in the polymer solution. This is because the density of segments near 
a point which is free of polymer is lower than in the bulk. The fractional reduction in density near 



a sphere/cavity is given by the term in parentheses on the right of Eq. ^ [20|. It is illustrated in 
Fig. 2. The methodology of the derivation of Eq. (^) is quite general, it may also be used to derive 
yp for ideal chains and for other dimensionalities. 

The positions of the spheres in the polymer solution are correlated over a range ^; they are 
distributed as if there is an attractive part to their interaction. The effective attraction is expressed 
as an effective potential between two spheres, called the potential of mean force w{r) [p. This is 
defined by 

w{r) = -T\ngs{r) (6) 
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which is, using Eq. (^), 



w{r) 



oo r < D 

-Tin 



1 I (^/O'^' .. (nn(r) 



r>D • 



As Z)/^ ^ 1, we may hnearise the logarithm 



w{r) ~ -T ^^^^ ^ (5p(r) - 1) r > i?. (8) 

(1 - 

Then from Eqs. (|l|) and (|8|) we see that the depth of the well in the potential of mean force is 
T{D / ^)^/^ <C T. The potential is weak but long ranged, and as the range is increased the depth of 
the well decreases. 

Now that the distribution function between spheres has been obtained we can calculate the 
structure factor 5(g), where q is the wavevector. S{q) could be measured in a scattering experiment 
performed on the mixture; note that it is not the total structure factor for the mixture, just the 
partial structure factor of the spheres. If we define the total correlation function for the spheres 
hs{r) = gs{r) - I [|ll],then 

S{q) = l + cshs{q). (9) 
The Fourier transform of the total correlation function, hs{q), is 

hs{q) = J hs{r)e'''-''dr, (10) 

where hp is the polymer segment-segment total correlation function. Using Eqs. (^) and (^), Eq. 
( |lO[ ) becomes 

. - f e'-* + (^^^g^ /; »Mr)e'^-'^r, (U) 

where hp{r) = gp{r) — 1. The first integral of Eq. ( |Tll) is just the lowest order, in density, 
approximation to the Fourier transform of the total correlation function between hard spheres 
The second term of Eq. ( pA] ) is not the Fourier transform of hp{r) because the region of integration 
is restricted to r > D. However, this region only contributes a fraction (D/^)^/^ to the integrand, 
for q <C D~^, and so we approximate the integral of the second term of Eq. ( [Tl| ) by hp{q). Then, 
Eq. (|Tll) becomes 

M.)-^-=(^).(:«^M., 

where ji is the first spherical Bessel function, ji[z) = sin(z)/2;^ — cos{z)/ z. We have related hs{q) 
and hence S{q) to the polymer segment-segment correlation function hp and so /ip, and therefore 
can be measured using only scattering from the spheres. 
For hp{q) we use the Ornstein-Zernike form, valid for small q, q ^ C^^, 

hpil) - (^^, (13) 

which is the Fourier transform of the third line of Eq. ([l|) . Although this gives the correct behaviour 
at small g, at g ^ ^, hp varies as not the q~'^ given by Eq. (13): 5/3 is the exponent for 



a self-avoiding walk and at length scales smaller than ^ the segment-segment interactions are not 
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screened and the chains are self-avoiding ||2|]. We take Eq. (|13D as unless the difference between ^ 
and D is very large, say O(IOO) which requires a very large polymer, then the behaviour at small 



enough values of q to distinguish between 5/3 and 2 will be obscured by first term in Eq. (12) and 
so not observable in measurements of S{q). Then, 



In Fig. 3 we have plotted hs{q) both with and without a semidilute polymer solution. The 
hsil) is just what we would expect for particles interacting via a short ranged repulsion — the 
hard sphere interaction — and a long ranged attraction — the attractive potential of mean force 
due to the polymer. Fig. 3 may be compared with the Fig. 5 of Khalatur et al. jl^ who observe 
the same effect. Note that at large values of q the polymer has no effect, hs is dominated at large 
q by the hard-sphere interactions. While at low q hs is less negative when polymer is present. We 
can compare this with a van der Waals fiuid in which the attractive interaction is infinitely long 
ranged; in that case hs{q) = hnsiQ), for hns the hard-sphere total correlation function, except for 



q = where hs > h^s by an amount proportional to the extent of the attraction [18]. Due to the 
long range of the attraction, S,, as compared to the repulsion, D, the spheres behave very much 
like a van der Waals fiuid. Indeed, they are an even better approximation to a van der Waals fiuid 
than are simple fiuids such as the noble gases. The ratio of the range of the attractive interaction 
to that of the repulsive interaction is smaller for noble gases than for the spheres in the polymer 
solution. However, this is only true at low densities of spheres. When the density of spheres becomes 
comparable to the density of polymer (rescaled) polymer segments the correlations between polymer 
segments will be affected by the spheres and E, reduced. The distribution of spheres in the mixture 
is much more uniform than that of the polymer, although as we have shown there are significant 
correlations in the positions of the sphere. So, we estimate that even if the the density of spheres 
is comparable to that of the polymer segments the spheres will not greatly affect the correlations. 
The situation where the density of spheres is high is quite subtle and so far beyond this work. 

The interactions studied here between spheres in a polymer solution are similar to those between 
spheres in a near-critical binary mixture |^]. There too the spheres behave as if they attracted 
each other, due to correlations (inhomogeneities) in the solvent in which they are suspended. Both 



these interactions are analogous to the quantum-mechanical Casimir interaction |21]. 

3 Conclusion 

The pair distribution function of small hard spheres in a semidilute polymer has been obtained, Eq. 
(|l4|), and it has been shown how to determine the polymer's correlation length from the structure 
factor of the spheres. The correlations in distribution function function are long ranged, i.e., the 
ratio of the range of the polymer induced effective attraction, ^, to the range of the repulsion, D, 
is much larger than the ratio of the attractive to the repulsive forces in simple fiuids such as the 
noble gases. Note that this effective attraction is in a mixture where all interactions are purely 
repulsive excluded volume interactions. As the range of the attraction is just ^, it is a function of 
the polymer concentration and can be varied at will within the semidilute regime. Also the depth 
of the attractive well between two spheres varies with ^; as the range increases the well depth 
decreases. Finally, we note that as only the ratio in size between the particle and the polymer's 
correlation length is important, the behaviour described above would also occur if the polymer is 



not a fiexible polymer but a semifiexible biopolymer such as DNA or wormlike micelles [p2| , 23 1 and 
the particles are much larger particles with diameters of a few hundred nm. 
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Figure Captions 



Fig. 1 A schematic picture of a small colloidal sphere, a micelle or a protein molecule, and a 
nearby part of a polymer chain. The black circle is a spherical particle. The chain is shown as a 
black curve and the dashed circles represent the rescaled segments of size D. 

Fig. 2 The polymer density around a sphere/cavity in the polymer solution. If there is a sphere 
at the origin then there can be no polymer segment there and so the polymer density is reduced in 
the surrounding volume. The x-axis of the graph is r the distance from the sphere and the y-axis is 
the polymer density. The segment density in contact with the sphere is less than cp but non-zero; 
it relaxes to the bulk value cp at a rate ~ r~^/^ for r < ^ and exponentially when r > ^. 

Fig. 3 The total correlation function for low-density hard spheres. The solid curve is for pure hard 
spheres and the dashed curve is for spheres in a semidilute polymer solution with ^/D = 10. The 
inset shows the small q part at a larger scale. 
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Fig.l 




Fig. 2 



